Critical behavior of YBa2Cu307 ^ and Bi2Sr2CaCu208+<5 from the dual 

Ginzburg-Landau model 
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An unifying scenario for the critical behavior of YBa2Cu307_5 (YBCO) and Bi2Sr2CaCu208+4 
(BSCCO) is proposed. It is shown that both critical crossovers observed in these materials follow 
by considering two different scalings in the dual Ginzburg-Landau model. The first scaling leads to 
the critical exponents v « 2/3, u' ~ 1/3 and a f» 0, with u and a being respectively the correlation 
lenght and specific heat exponents while z^' is the magnetic field penetration depth exponent. These 
values for the critical exponents agree with the ones obtained experimentally for YBCO single 
crystals. For the second scaling we obtain i/ = 1 and a = —1 which must be compared with the 
measured values f Ri 1 and a « —0.7 for BSCCO. For the penetration depth exponent it is predicted 
the value i/' = 1. 
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It is generally accepted that the experimentally acces- 
sible critical region of high temperature superconductor 
YBa2Cu307_5 (YBCO) lies in a ^He universality class 
This means that the critical behavior is governed 
by the 3D XY {XY^ for short) nontrivial fixed point. The 
XY3 critical behavior of YBCO has beenprobed both in 
zero and nonzero external field regimes |3|. For instance, 
a quantity which can be accurately measured in zero field 
is the magnetic field penetration depth. A, whose scaling 
behavior near Tc is given by A \t\''' , with v' = 0.33±0.1 
lH and t being the reduced temperature. This value of 
ly' is consistent with a critical regime described by a XYt, 
universality class where we expect = 1^/2 with v w 2/3 
|4|. Also, specific heat measurements give a « —0.013 
which is also consistent with the XY3 universality 
class. These results hold for single crystals of optimally 
doped YBCO. However, this scenario may not hold for 
thin films of YBCO. For example, the value i'' — 1/2 has 
been measured at zero field regime which may be 
associated either to a mean-field behavior |^ or to some 
other critical crossover (7||^. It is worth to cite also the 
recent results of Charalambous et al. ^ for bulk YBCO. 
They reported different values of a depending on whether 
Tc is approached from below or from above. 

The XY3 universality class is not always shared by 
other cuprate superconductors. For instance, a dif- 
ferent critical behavior is probed in Bi2Sr2CaCu208+a 
(BSCCO). It is obtained that a « -0.7 and « 1 
This is not very different from the critical exponents of 
the 3D 0{N) model in the large N limit where we find 
exactly a = —1 and v = \ ||l2|. In the context of super- 
conductors, these exponents are obtained in a Hartree 
approximation in external field but with the gauge field 
fiuctuations being neglected Indeed, the Hartree ap- 
proximation consists of a gap equation analogous to the 
one obtained in the large N limit of the 0{N) model. 
The critical exponents a = —1 and j/ = 1 can be ob- 
tained also through a weakly interacting Bose gas model 
|[14[. In this case the specific heat exhibits a triangular 



peak at zero field characteristic of a Bose-Einstein con- 
densation (BEC). 

In this paper we will show that the critical crossovers 
observed in both YBCO and BSCCO are described by 
two different scalings in the dual version of the Ginzburg- 
Landau model ^,|l^,^ (to be referred from now on as 
the dGL model). The description of the XY^ scaling of 
YBCO by the dGL model is well known jl|,||. How- 
ever, we want to emphasize the fact that two different 
crossovers are obtained here from a same model. We 
provide in this way an unifying view for the YBCO and 
BSCCO critical behaviors. 

The dGL model was introduced in the literature as a 
continuum version of the lattice dual Ginzburg-Landau 
(GL) model ||l^. The dGL model is therefore a way to- 
wards a field theoretical description of the inverted XY3 
transition ||l^ . The inverted XY3 scenario must be valid 
at least for superconductors in the type II regime and en- 
sures the existence of a charged infrared stable fixed point 
for the GL model |l9|-p^. One important difference of 
the inverted XY3 relative to the ordinary XYt, universal- 
ity class is the scaling A ^ ^ m2l|,|2^ , where ^ is the cor- 
relation lenght. This scaling implies that u' = v instead 
of v' = v/2 as in the ordinary XY^ fixed point. Recent 
high precision Montecarlo simulations give further sup- 
port to the duality scenario and to the existence of the 
inverted XY^ critical point j2^-|2^]. Unfortunately, the 
critical region corresponding to this nontrivial charged 
fixed point is still experimentally out of reach. 

The bare free energy density for the dGL model is given 

by 



^=i[(Vxho)2 + A/oX] 



(V - ieo,dho)V'o| 



+ yl^ol 



(1) 



where the constraint V • ho = should be understood. 
The bare dual charge eQ,d = 27r A/0/90 where qo = 2eo is 
the Cooper pair charge. Also, — 9o(l'/'oP( where (/)o 
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is the bare order parameter field of the GL model. Thus, 
Mq is the photon mass generated in the Meissner phase 
of the GL model. The Meissner phase of the GL model 
is described as a symmetric regime in the dGL model. In 
this sense, the dGL model is a disorder field theory. 

Due to the presence of two masses in the problem, we 
have more possibilities of scalings in the dGL model than 
in the GL model. A class of such scalings was discussed 
recently by us Q . 

The scaling corresponding to the XY^ behavior of 
YBCO is obtained by looking to the fiow of the renor- 
malized dimensionless couplings / ~ e^/m and g = u/m 
with respect to m at Mq, uq and eo,d fixed. The renormal- 
ized parameters are defined in a usual way . They are 
given by = Z-^Z^ml, = ZhM§, = Zhcl ^ and 
u = Z'^uq/Zu- Such a scaling has been considered before 
in Refs. and for this reason we will omit the de- 

tails. In this scaling the infrared stable fixed point corre- 
sponds to a nonzero /. As a consequence, the anomalous 
dimension rjfi = mdlnZii/dm is given at the nontrivial 
fixed point by rj^ = 1 exactly Since Mq is fixed 

it follows the scaling M ~ m^'^near the critical point. 
Thus, it follows that ly' — vjl and the ratio m/M — > 
as m ^ 0. From this last observation we obtain that the 
gauge degrees of freedom decouple and the correspond- 
ing fixed point is XY^. We have therefore that v w 2/3, 
v' « 1/3 and a « 0, which are just the measured values 
in the bulk YBCO with optimal doping 

For the scaling corresponding to the BSCCO critical 
behavior it is more convenient to choose M as the run- 
ning scaling variable. The BSCCO scaling is defined by 
demanding that the ratio jtiq/Mq is kept fixed together 
with eo,(i and uq We assume therefore that both 

mp and Mq are proportional to the reduced tempera- 
ture t. In this scaling the renormalized quantities are 
defined as before and the corresponding rcnormalization 
constants are denoted by Z^, Z^^ Z,n and Z„. The di- 
mensionless couplings are now defined by / = ZhS^ JM 
and g = Z~^Z|,Mo/Af. It is convenient to define also the 
dimensionless parameter = m/M which plays a role 
analogous to the Ginzburg parameter. Note that fixing 
eo,d is equivalent to consider that the amplitude fiuctua- 
tions are frozen in the original GL model. We define the 
rcnormalization group (RG) functions: 



M 



d\nZ, 
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dM ' 
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We have the following exact fiow equations: 



(2) 
(3) 
(4) 

(5) 



M 



dm? 
'dM 



(6) 



(7) 



(8) 



A nontrivial fixed point with /, g and nonzero must 
verify the equations fjh — 1, fj^ — fjm — 

1 and Mdg/dM = 

0. If such a fixed point exists, we obtain from Eqs. (Q) 
and (^) that v' = v = 1 and a = —1 from the scaling 
relation 31^ = 2 — a, which are the required exponents. 
Note that by assuming the existence of the nontrivial 
fixed point we have that the values of the exponents are 
exact and are the same as in the BEG transition. The 
existence of this nontrivial fixed point can be verified ap- 
proximately through a simple 1-loop example. We have 
at 1-loop order. 
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where in neglecting higher order terms we have assumed 
Kd,o — ttio/Mq « 1. This assumption does not affect the 
generality of the problem though higher values of Kd^ 
lead to a worse result from the point of view of conver- 
gence. The numerical nontrivial fixed point values are 
« 0.34 and g* ^ 21.4. The numerical fiow diagram 
in the gK^-plane is shown in Fig. 1. This result can be 
improved further by computing higher order terms and 
using resummation methods. Although we have shown 
the validity of the proposed scaling through perturbative 
methods, we claim that this scenario holds also beyond 
perturbation theory. As mentioned earlier, the present 
scaling with cq and mo/Mo fixed corresponds to freeze 
the amplitude fiuctuations in the original GL model. Nu- 
merical simulations with frozen amplitudes but taking 
into account the phase fluctuations can be used as a final 
test of the above picture. 

Note that the scaling we have just discussed predict 
the value v' = I for the penetration depth exponent. 
This prediction implies the scaling A ~ ^ which is veri- 
fied also in the vicinity of the inverted XY^ critical point. 
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It would be very interesting to verify this prediction ex- 
perimentally and also by numerical means. 

Summarizing, we have considered two different scalings 
in the dGL model, each one corresponding to a different 
critical crossover. The first one gives the XY3 scaling be- 
havior which is probed experimentaly in single crystals 
of YBCO while the second one reproduces approximately 
the scaling behavior of BSCCO. Moreover, we predict 
the value ly' — 1 for the BSCCO penetration depth ex- 
ponent. The results of this paper show the power of the 
dual approach. Indeed, all the possible crossovers ob- 
tained within the dGL model may be very difficult to 
obtain in the original GL model. The point is that a 
weak coupling regime in the dGL model corresponds to 
a strong coupling regime in the GL model. In this case 
nonperturbative effects become important. 

Among the further directions is the inclusion of 
anisotropy and more general disorder parameters with 
higher symmetries like d + is wave [|9| and 50(5) ||3C| ] 
symmetries. It would be also interesting to study other 
nouvel situations which occur in external field like the 
^-transition of Tesanovic . Finally, we hope that the 
present contribution will stimulate some experimental ef- 
fort towards the verification of the prediction v' = 1. 

We would like to thank Z. Tesanovic for bringing Ref. 
1^ to our attention and A. Sudb0 for sending us Ref. |2^] 
prior publication. We would like to thank specially Prof. 
A. Junod for his valuable comments on the experimental 
results in YBCO and BSCCO. 
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